home *** CD-ROM | disk | FTP | other *** search
/ Celestin Apprentice 4 / Apprentice-Release4.iso / Source Code / Libraries / VideoToolbox 95.11.08 / Utilities / Quick3 / TestPsychometricFit.c < prev    next >
Encoding:
Text File  |  1995-06-17  |  4.3 KB  |  33 lines  |  [TEXT/CWIE]
  1. m.\n",(int)weibullDF);
  2.     printf("Observer: ");
  3.     printf("logAlpha%6.2f, beta%4.1f, gamma%5.2f, delta%5.2f\n",params.logAlpha,params.beta,params.gamma,params.delta);
  4.     p=PsychometricFit(¶ms,&Weibull,&data,&weibullLL,weibullDF,&chiSquare,&chiSquareDF);
  5.     printf("Fit:      ");
  6.     printf("logAlpha%6.2f, beta%4.1f, gamma%5.2f, delta%5.2f, significance%5.2f\n",params.logAlpha,params.beta,params.gamma,params.delta,p);
  7.     #if QUEST
  8.         printf("QUEST mode %.2f\n",mode);
  9.     #endif
  10.  
  11.     /*
  12.     We're done, but just to show off, let's print out everything that anyone
  13.     could possibly want. In real life I would skip this junk.
  14.     */
  15.     monotonicData=data;
  16.     MonotonicFit(&monotonicData,&monotonicLL,&monotonicDF);    /* overwrites data with fit */
  17.     printf("\ncontrast Trials Right  Ratio  Weibull Monotone\n");
  18.     for(i=0;i<data.contrasts;i++){
  19.         cPtr=&data.c[i];
  20.         printf("%6.3f   %5ld %5ld %7.3f %7.3f %7.3f\n",
  21.             cPtr->contrast,cPtr->trials,cPtr->correct,
  22.             cPtr->correct/(double)cPtr->trials,
  23.             Weibull(cPtr->contrast,¶ms),
  24.             monotonicData.c[i].correct/(double)monotonicData.c[i].trials
  25.             );
  26.     }
  27.  
  28.     chiSquare=2.0*(monotonicLL-weibullLL);
  29.     chiSquareDF=monotonicDF-weibullDF;
  30.     p=PChiSquare(chiSquare,chiSquareDF);
  31.     printf("\nChi square %.1f with %d degrees of freedom, yielding a significance of %.2f\n"
  32.         ,chiSquare,(int)chiSquareDF,p);
  33. }